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I.  INTRODUCTION 


To  properly  describe  coupled  chemical 
reactions  and  gaseous  diffusion  in  porous 
sorbent  and  catalyst  grains,  the  "pore  tree"  was 
introduced  by  Simons  and  Finson  (1979)  and 
Simons  (1982).  The  pore  tree  represents  an 
isolated  sub-structure,  allowing  diffusion  into 
and  out  of  the  porous  media  without  permitting 
transport  through  the  media.  This  f>ore  structure 
was  develop>ed  via  analogy  to  the  kinetic  theory 
of  gases  with  the  pore  length  analogous  to  the 
mean  free  path.  Under  the  assumption  that  the 
pore  aspect  ratio  (pore  length  to  radius)  is  a 
constant,  a  pore  size  distribution  was  obtained 
that  has  been  confirmed  for  coal,  coal  char, 
sorbents,  catalysts  and  kidney  stones  from  both 
men  and  women.  The  pore  tree  was  statistically 
derived  from  the  pore  size  distribution  and 
allows  the  orderly  migration  of  a  reactant  gas 


Figure  1 .  The  Pore  Tree 


from  the  large  pores  to  the  small  pores  (Fig.  1). 

A  detailed  description  of  the  pore  tree  and  the  coupled  transport  and  chemistry  is  given  by 
Simons  (1982,  1983a).  The  spatially  dependent  transport/reaction  equations  are  solved  for  a 
single  pore  tree  and  then  the  total  contribution  of  all  trees  (of  all  sizes)  in  the  system  is 
obtained  by  summing  the  contribution  of  each  tree  that  reaches  the  exterior  of  the  system. 

This  is  distinct  from  the  "bulk"  transport  approach  in  which  the  transport  equation  for  a  single 
pore  is  integrated  over  all  pores  at  a  fixed  point  in  space  before  integrating  spatially.  The 
"bulk"  transport  approach  is  invalid  if  the  spatial  gradients  in  the  transport  equations  are 
implicit  functions  of  pore  size.  One  example  of  this  implicit  pore  size  dependence  is  that  of 
the  heterogeneous  reactions  within  porous  catalysts  and  sorbents  for  which  the  pore  tree 
structure/transport  model  was  developed.  A  second  example  is  that  of  coupled  diffusion  and 
remediation  reactions  in  the  immobile  region  of  soil. 


In  order  to  describe  the  subsurface  transport  of  gas  and  water  in  soil,  the  dispersion  of 
contaminants,  and  in-situ  remediation  of  contaminated  sites,  the  pore  tree  is  extended  herein 
to  simulate  permeability  and  bulk  transport.  The  random  nature  of  the  pore  structure,  which 
formed  the  basis  of  the  statistical  derivation  of  the  pore  tree,  is  applied  to  porous  soil  and 
sand.  The  interconnectivity  of  the  pore  structure  is  obtained  via  a  statistical  determination  of 
the  "branches"  that  are  common  to  several  trees  to  allow  convection  and  diffusion  through  the 
large  scale  (mobile)  structure  in  addition  to  diffusion  and  coupled  reactions  in  the  small  scale 
(immobile)  structure.  The  extended  pore  tree  model  has  been  used  to  explain  the  measurement 
errors  in  the  permeability  of  soil  due  to  the  measurement  scale  size  (Shouse  et.  al.,  1994)  and 
has  successfully  predicted  the  bulk  gaseous  diffiisivity  m  partially  saturated  soil  (Washington 
et.  al.,  1994)  as  a  function  of  a  saturation  scale  size.  The  theory  also  depicts  a  permeable 
sub-range  in  which  small  scale  convection  bridges  the  flow  between  the  bulk  convection  and 
the  small  scale  diffusion.  The  model  provides  an  analytic  description  of  a  pore  structure  for 
soil  upon  which  transport  and  coupled  chemical  reactions  may  be  accurately  superimposed. 
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n.  ISOLATED  PORE  TREE:  THE  STRUCTURE 


Following  the  pore  structure  theory  of  Simons  and  Finson  (1979)  and  Simons  (1982), 
consider  a  spherical  porous  particle  of  radius  a,  containing  pores  of  length  Ip  and  radius  rp. 
The  pore  dimensions  range  from  a  microscale  of  the  order  of  Angstroms  to  a  macroscale 
which  is  a  significant  fraction  of  the  particle  radius.  The  radius  of  the  largest  pore  is  denoted 
by  r„,„  and  is  given  by 


r^-2ae'P/3K^  (» 

where  6  is  the  total  porosity  of  the  particle  and  is  a  constant  of  integration, 
approximately  equal  to  five,  which  relates  the  pore  length  to  its  radius 

The  radius  of  the  smallest  pore  is  denoted  by  r„,in  and  is  given  by 


'•.*.=2e/Pp.s,  (3) 

where  ps  is  the  density  of  the  solid  matrix,  Sp  is  the  specific  internal  surface  area  (several 
hundred  mVg),  and 


The  particle  contains  a  continuous  distribution  of  pore  sizes  from  r„,i„  to  r,„j„  •  The 
number  of  pores  within  an  arbitrary  plane  of  cross-sectional  area  A  and  with  radius  between 
rp  and  rp  +drp  is  denoted  by  g(r^)Adr^.  The  pore  distribution  function  5(r^)is  given  by 


8irj,)=Q/2%^r^  (5) 

where  g(r^)  indicates  an  average  over  all  inclination  angles  between  the  axis  of  the  pore  and 
the  normal  to  the  plane.  Due  to  the  random  orientation  of  the  pores,  the  intersection  of  a 
circular  cylinder  with  a  plane  is  an  ellipse  of  average  area  27trp^.  Hence,  the  porosity  is  the 
27t:rp^  moment  of  ^(r^)  and  the  internal  surface  area  is  the  47irp  moment  of  ^(r^) .  The 
expression  for  gir^)  was  derived  (Simons  and  Finson,  1979)  from  statistical  arguments  and 
has  been  validated  through  extensive  comparison  of  the  predicted  volume  and  surface  area 
distributions  with  mercury  intrusion  data  (Stacy  and  Walker,  1972).  This  has  been 
accomplished  for  coal,  char  derived  from  that  coal  (Kothandaraman  et.al.,  1984),  sorbents, 
catalysts  and  even  kidney  stones  from  both  men  (low  porosity  oxalate)  and  women  (high 
porosity  phosphate). 
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A  characteristic  feature  of  the  distribution  depicts  that  the  pore  volume  between 
r^in  and  tp  increases  linearly  with  the  natural  log  of  rp.  It  is  the  functional  form  of  this 
relationship, 


Pore  Volume  “  J Si^p)drp  «lnr^ 


uiin 


(6) 


that  depicts  the  inverse  cubic  dependence  of  ^(r^)  on  rp.  A  linear  display  of  mercury 
intrusion  volume  vs.  In(rp)  always  infers  a  Mx^  distribution. 

The  number  of  pores  within  the  bulk  volume  V  whose  pore  radius  is  between  x^  and 
tp+dtp  may  be  defined  by  Vfir^)dr^.  The  pore  volume  is  expressed  as  the  tttp^lp  moment  of 
/(r^)  and  the  internal  surface  area  is  the  27crplp  moment  of  f(r^) .  The  pore  size  distribution 
functions  ( /(r^)  and  g(r^)  )  are  clearly  not  independent.  The  definitions  of  porosity  and 
internal  surface  area  infer  that  f(r^)  is  related  to  ^(r^)  by 


g(rp)  =fir^)lpl2 


(7) 


Equation  (7)  simply  states  that  the  probable  number  of  pores  intersecting  an  arbitrary  plane 
increases  with  the  length  of  the  pore  and  with  the  density  of  pores. 


The  length  of  a  pore  is  determined  by  an  arbitrary  intersection  with  another  pore  and 
is  expressed  (Simons  and  Finson,  1979)  as  a  collision  integral  over  the  pore  distribution 
functions.  The  analysis  suggests  that  Ip  ,  g'(r^)  and  /(r^)  are  proportional  to  rp ,  lAp^  and 
l/Xp  respectively.  The  constants  of  proportionality  are  obtained  from  integral  constraints,  i.e., 
the  total  porosity  and  internal  surface  area  contained  in  the  pore  structure.  The  expression  for 
is  given  by 


f(rp)  = 


(8) 


where  the  constants  were  defined  above. 

The  pore  volume  distribution  corresponding  to  these  distribution  functions  is  similar  to 
that  utilized  in  the  random  pore  model  (Gavalas,  1980  &  1981).  However,  the  pore  tree 
model  and  the  random  pore  model  differ  dramatically  in  their  choice  of  the  pore  aspect  ratio 
(length  to  diameter)  and  its  implications  with  respect  to  pore  branching.  The  random  pore 
model  allows  a  single  pore  to  connect  two  larger  pores.  This  picture  lends  itself  to  the 
idealization  of  instantaneous  mixing  between  the  pores  and  requires  that  the  pore  aspect  ratio 
be  of  the  order  of  one  hundred.  The  pore  tree  theory  uses  data  for  to  imply  (via  KJ  that 
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all  pores  possess  an  aspect  ratio  of  the  order  of  ten.  Hence,  small  pores  may  connect  to  larger 
pores  only  on  one  end  and  all  pores  must  branch  from  successively  larger  pores  like  a  tree  or 
river  system. 

Each  pore  that  reaches  the  exterior  surface  of  the  particle  is  depicted  as  the  trunk  of  a 
tree.  The  size  distribution  of  tree  trunks  on  the  exterior  surface  of  the  particle  is  denoted  by 

g{r^)A%a^dr^  where  g{r^  is  functionally  identical  to  g(r^).  Each  trunk  of  radius  r,  is 
associated  with  a  specific  tree-like  structure.  Let  N,  be  defined  as  the  branch  distribution 
function  where  N,drp  is  the  number  of  pores  of  radius  rp  (within  size  range  dtp)  in  a  tree 
whose  trunk  radius  is  Tj  .  The  total  number  of  pores  of  radius  tp  in  a  sphere  of  radius  a 

may  be  expressed  as  4l3na^f(r^)dr^  or,  as  the  sum  of  all  pores  of  radius  rp  contained 
within  every  tree  in  the  porous  sample,  plus  all  pores  of  radius  rp  that  are  themselves  the 
trunk  of  a  tree.  Hence, 


=  J  N^g(r^)4v:a^dr^  +  4na^g(r^)  (9) 


where  g(r^)  is  the  number  of  tree  trunks  per  unit  external  area  of  the  porous  sample  and  only 
those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp .  Using  the 
previously  derived  expressions  for  r„,„,  g(r^)and  f(r^) ,  Eq.(9)  is  identically  satisfied  by 


XT  3,4 


(10) 


The  branch  distribution  function  completely  characterizes  the  pore  tree.  The  internal 
surface  area  and  pore  volume  associated  with  each  pore  tree  are  denoted  by  St(r,)  and  V,(rJ, 
respectively,  and  are  expressed  as  the  sum  of  the  contributions  from  the  trunk  and  that  from 
the  branches. 


5,(0  =2itr,/,  +  j  2%r^l^N,drp 


(11) 


1^,(0  =nr^l^+fi:r^ 


(12) 


*  mtn 


Using  Eq.(lO)  for  N,,  St(rJ  and  Vt(ri)  become 
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(13) 


(1-6) 


/ 

l+ln 

t 

\ 

/min)) 

(14) 


where  the  (1-0)  term  in  St  has  been  included  to  account  for  pore  combination  (Simons,  1979). 

The  surface  area  associated  with  the  pore  tree  may  be  several  orders  of  magnitude 
greater  than  the  surface  area  of  the  trunk.  However,  the  volume  of  the  pore  tree  may,  at 
most,  be  one  order  of  magnitude  greater  than  that  of  the  trunk.  It  should  also  be  noted  that 
the  above  expressions  for  S^  and  Vj  reduce  to  those  appropriate  to  a  single  cylindrical  pore  in 
the  limit  of  rt  — >r^  (the  leaf  of  the  tree).  Furthermore,  the  integrals  of  St(rt)  and  V/rJ  over 
all  g(rj)  recover  the  total  internal  surface  area  and  pore  volume  of  the  porous  sample. 

Each  trunk  of  radius  rj  is  associated  with  a  specific  tree-like  structure  with  continuous 
branching  to  ever  decreasing  pore  radii.  The  radius  and  number  of  pores  is  a  unique  function 
of  the  distance  x  into  the  tree.  The  coordinate  x  is  skewed  in  that  it  follows  a  tortuous  path 
through  the  branches  of  the  tree.  Let  n(x)  represent  the  number  of  pores  of  radius  rp  at 
location  x  in  a  tree  of  trunk  radius  q.  An  analysis  (Simons,  1982)  of  this  pore  tree  has 
demonstrated  that 


n(x)=rflrj(x) 


(15) 


and  the  coordinate  x  is  related  to  rp  by 


dr^ldx  =  -r^ll^ 


(16) 


The  continuous  branching  model  has  been  used  to  successfully  describe  char  oxidation 
(Lewis  and  Simons,  1979;  Simons,  1982  &  1983a),  coal  pyrolysis  (Simons,  1983b  &  1984) 
and  the  catalytic  cracking  of  benzene  by  porous  iron  oxides  (Simons  et  al.,  1986).  It  was  also 
used  to  successfully  describe  sulfur  sorption  (SOj  and  HjS)  by  porous  calcine  (CaO)  in  the 
limit  of  zero  utilization  (Simons  and  Rawlins,  1980;  Simons  et  al.,  1984)  and  was  later 
extended  to  include  CaSO^  and  CaS  deposits  (Simons  and  Carman,  1986;  Simons  et  al.,  1987; 
Simons,  1988;  Simons  et  al.,  1988).  The  subsequent  determination  of  the  controlling  physical 
parameters  led  to  a  new  concept  for  the  optimization  of  the  sulfur  sorption  process  (Simons, 
1991;  Simons  et  al.,  1992)  through  spray  drying  of  water  soluble  organic  calcium  solutions  to 
control  the  sorbent  pore  structure. 
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in.  INTERCONNECTIVITY 

The  first  step  in  determining  the  size  distribution  of  the  interconnected  pores  and  the 
distribution  of  the  permeability  is  to  determine  the  distribution  function  Gf(r^,r^dr^  which 
represents  the  number  of  pores  of  radius  rp  (within  size  range  drp)  per  unit  cross  section  of  an 
arbitrary  plane  and  also  contained  within  a  tree  whose  trunk  radius  is  r^  Consider  an  infinite 
homogeneous  isotropic  porous  media  and  isolate  a  spherical  volume  of  that  media  denoted  by 
the  radius  a.  Such  a  volume  is  illustrated  in  Fig.  2.  The  total  number  of  pores  of  radius  rp 
(within  size  range  drp)  intersecting  plane  AA  of  area  Jia^  has  previously  been  defined  by 

g(r^)na^dr^  .  The  pores  in  plane  AA  in  this  size  range  may  also  be  determined  by 

integrating  Gf(r^,r^T:a^dr^  over  all  trees  whose  trunk  intersects  the  exterior  surface  of  the 
porous  sample.  Hence  it  follows  that 


g(r^) It a^dr^  =  J  [G,(r^,rp n a^dr^]  g(r^) 4na^ dr^  (1'7) 

0- 

where  only  those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp . 

A  solution  to  Eq.  (17)  for  G,(rj,rp  will  not  necessarily  be  unique.  Physical  arguments 

will  help  determine  G,(r^,rp  and  help  ensure  that  it  is  the  particular  solution  we  seek.  Since 
N,  represents  the  number  of  pores  of  size  rp  in  the  tree  and  the  probability  of  a  pore 

intersecting  a  plane  is  proportional  to  its  length,  it  follows  that  G,(r^,rp  should  be 


Figure  2.  Spherical  Volume  of  a  Porous  Media 
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proportional  to  the  product  of  Nt  and  Ip/l^ ,  i.e.,  proportional  to  Eq.  (17)  is  identically 

satisfied  by  a  function  which  differs  from  r^k^  by  In(rp). 


4:ra^r^ln(r 


max 


/r,) 


(18) 


Note  that  ln(rp/r^)  introduces  an  integrable  singularity  at  rp=r^  such  that  0^{r^,r^)dr^  is 
finite  at  rp=r^  .  Hence,  there  is  one  and  only  one  largest  pore  for  each  reference  sphere. 

The  probability  of  trees  sharing  common  branches,  i.e.,  the  interconnectivity  of  the 

pore  structure  is  described  in  Fig.  3.  We  seek  the  distribution  function  I{rp)dr^  which 
represents  the  number  of  pores  (within  size  range  drp  about  rp)  per  unit  area  of  plane  AA  that 
are  connected  to  both  sides  of  the  pore  structure  through  pores  at  least  as  large  as  rp.  A^  is 
defined  as  the  area  within  plane  AA  that  is  open  to  one  side  of  the  porous  media  through  all 
trees  of  size  (through  all  pores  of  size  rp’  that  are  at  least  as  large  as  rp).  Subsequently, 

Ag  Gf(T^,T^drp  represents  the  number  of  pores  of  size  rp  (within  size  range  drp)  per  unit  area 
of  plane  AA  that  are  contained  in  a  tree  of  size  range  drj  about  r,  and  are  also  connected  to 
the  opposite  side  of  the  porous  media  through  all  trees  denoted  by  r,' .  It  follows  that  the 
distribution  function  for  interconnected  pores  in  plane  AA  may  be  obtained  by  integrating 

Ag  Gf{r^,r^dr^  over  all  trees  (r,)  that  are  large  enough  to  contain  a  pore  of  size  rp  .  Hence, 


I{rp)%a^dr^  =  f  lAQG^(r^,r^dr^  gir^)2Tza^dr^ 


(19) 


Figure  3.  Interconnectivity  of  a  Porous  Media 
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From  the  above  definition  of  may  be  expressed  as 


^6  ■  /[/ ^o^G,ir„r^dr^  g(r^)2%a^dr, 


(20) 


where  the  primes  on  the  variables  of  integration  have  been  omitted.  Evaluating  Eq.  (20) 
yields 


2p 


from  which  Eq.  (19)  yields  the  common  branch  distribution  function. 


(21) 


I(r^  = 


0  /r  ) 


4p 


(22) 


It  has  been  deduced  that  the  total  number  of  common  branches  of  size  rp  in  an 
arbitrary  plane  scales  approximately  with  the  total  number  of  pores  of  that  size  in  that  plane. 
Hence,  there  is  a  probability  of  interconnectivity  at  pore  size  rp  that  is  logarithmic  in  pore 
size.  Defining  this  probability  as  Pi(rp)  via  Eq.  (22), 


4p 


(23) 


it  is  apparent  that  approximately  one  percent  of  all  pores  of  all  sizes  are  interconnected 
through  larger  pores. 

The  broad  size  range  associated  with  the  interconnectivity  suggests  that  a  very  wide 
range  of  pore  sizes  control  transport  and  that  a  complicated  mixture  of  convective  and 
diffusive  transport  persists  through  all  of  pore  space.  While  permeability  is  dominated  by  the 
largest  pores,  it  is  important  to  determine  the  level  of  convection  that  is  occurring  in  smaller 
pores  in  order  to  accurately  describe  the  fine  scale  transport  necessary  to  assess  chemical 
reactions.  In  the  following  sections,  bulk  permeability,  bulk  diffusivity  and  small  scale 
convection  are  addressed  using  the  interconnectivity  derived  above. 
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IV.  PERMEABILITY 


The  fundamental  relationship  governing  convection  in  a  porous  media  is  Darcy's  Law 


which  relates  the  volume  flow  rate  to  the  pore  radius  rp  and  the  pressure  gradient 


dx 


(24) 


where  |i  is  the  viscosity  of  the  fluid.  The  bulk  permeability  (k)  is  defined  in 
terms  of  the  volume  flow  rate  across  the  cross  sectional  area  A 


A{-dpldx) 


(25) 


Convection  across  plane  AA  in  Fig.  4  will  possess  contributions  from  two  primary 
sources  illustrated  in  Figs.  4a  and  4b.  Fig.  4a  illustrates  the  case  where  the  convection  in 
plane  AA  is  due  solely  to  the  pores  that  are  interconnected  in  that  plane.  Fig.  4b  illustrates 
the  case  where  the  convection  in  plane  AA  is  due  to  the  smaller  pores  in  the  pore  tree  that 
are  interconnected  outside  of  plane  AA.  This  connectivity  will  translate  into  a  slower  velocity 
in  the  pore  crossing  plane  AA  but  could  be  significant  because  99%  of  the  pores  in  plane  AA 
are  not  interconnected  in  that  plane. 


a)  Pores  Interconnected  in  Plane  AA  b)  Pores  Interconnected  Out  of  Plane  AA 

Fig.  4  Convection  in  the  Pore  Tree 
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Consider  any  pore  of  radius  r,  in  plane  A  A  of  Fig.  4b  to  be  the  trunk  of  a  tree.  Each 
pore  of  size  rp  within  the  tree  possesses  the  probability  Pjir^)  of  being  interconnected  and 

each  interconnected  pore  in  the  tree  will  carry  volume  flow  rate  Qp(fp)-  Since  there  are 
Ngdrp  (Eq.  10:  Ns=r3Vrp‘‘)  pores  in  size  range  drp  within  the  tree,  the  total  volume  flow  rate 

through  trunk  r^  in  plane  A  A  becomes 

^3 

QJO  =  /  <?,(<>)  N,  dr^  (26) 


or,  to  first  order, 


6ttr;^ln(r^/rJ/  dp\ 

32  pp  \  dx) 


+  H,ojr 


(27) 


Within  this  approximation,  it  is  seen  that  is  identical  to  the  volume  flowing 

through  the  pores  that  are  interconnected  within  plane  AA.  i.e.. 


(28) 


which  demonstrates  that  all  volume  flow  through  plane  AA  in  pore  size  r^  is  dominated  by  the 
interconnectivity  of  size  r,  in  plane  AA  and  not  by  the  interconnectivity  of  smaller  pores  in 
subsequent  branches  of  the  pore  tree.  Simply  stated:  case  4a  dominates  case  4b. 

Since  all  volume  flow  through  plane  AA  is  limited  by  the  interconnectivity  of  the 
pores  in  that  plane,  Eq.  (24)  for  the  volume  flow  rate  may  be  rewritten  to  include  all 
interconnected  pores  in  area  A.  Hence, 


(29) 


where  /(r^)  is  the  "common  branch  distribution  function"  given  by  Eq.  (22).  The  bulk 
permeability  (k)  is  then  expressed  as 


k  = 


It 

8 


(30) 
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upon  integration,  Eq.  (30)  becomes 


k  = 


0r„ 


\2 


16  p 


(31) 


Equation  (31)  resembles  a  dozen  other  expressions  available  in  the  literature  wherein  it 
is  concurred  that  the  bulk  permeability  is  dominated  by  the  largest  pores  in  the  media  but  the 
unknown  value  of  that  permeability  is  simply  replaced  by  an  unknown  pore  size  to  the  second 
power.  Since  the  pore  size  distribution  function  will  be  least  accurate  at  the  extreme  end  of 
the  size  range,  i.e.  at  r^^^x »  no  claim  can  possibly  be  made  that  the  numerical  constants  in 
Eq.  (31)  are  in  any  way  superior  to  those  derived  elsewhere.  One  important  advantage  of  the 
extended  pore  tree  model  is  that  it  characterizes  the  distribution  of  permeability  in  pore  space, 
a  feature  that  will  be  important  in  describing  fine  scale  contartiinant  transport  and  in-situ 
remediation.  A  second  advantage  is  the  ability  to  assess  statistical  errors  in  the  measurement 
of  the  permeability  as  a  function  of  the  measurement  scale  size.  This  exercise  is  also  a  good 
test  of  the  extended  pore  tree  model. 

Consider  a  soil  sample  with  the  following  physical  characteristics: 


Conductivity: 

Permeability: 

Porosity: 

Pore  Aspect  Ratio: 
ln(r^/r„iJ 

From  Eq.  (31),  it  follows  that 

Tmax  =  300  Fm 

and  subsequently  r„i„  =  20  A,  The  size  of  the  smallest  pore  is  not  an  important  parameter  for 
this  application  but  may  be  readily  adjusted  through  a  minor  variation  in  the  value  of  (3  (e.g., 
for  r,ni„  of  the  order  of  100  A,  p  =  10).  The  bold  assertion  made  in  applying  this  pore 
structure  model  to  soil  is  that  the  pore  size  distribution  is  valid  between  r^  and 

To  investigate  the  role  of  the  measurement  scale  size  on  permeability,  consider  the 
largest  pore  r„,„  contained  in  the  spherical  sample  of  radius  "a"  as  given  by  Eq.  (1).  Since 
r„„,  in  Eq.  (31)  for  the  permeability  represents  the  largest  pore  in  the  media,  the 
corresponding  value  of  "a"  is  denoted  a„,^  and  represents  the  largest  sample  size  for  which 
the  pore  sizes  will  scale  with  the  dimensions  of  the  sample.  From  Eqs.  (1)  and  (31) 


V  =2  cm/hr 
k  =0.6  Darcy 
0  =50% 

Ko=5 
P  =12 


a. 


max 


24g,p 

04/3 


(32) 


Each  sphere  of  radius  a„,a,t  will  contain  one  pore  of  size  r^^.  A  20  x  20  grid  of  these  spheres 
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will  be  characterized  by  the  dimension  40a^  and  contain  400  pores  of  size  r^.  Each  of 
these  pores  possess  probability  Pi(rp)  of  being  interconnected.  Following  Eq.  (23),  Pi(rp)  is 
approximately  0.0025  for  rp  sufficiently  close  to  Hence,  only  one  of  the  4()0  Ingest 
pores  in  this  20  x  20  grid  will  be  interconnected  and  the  error  in  the  measurement  of  the 
permeability  will  correspond  to  the  statistical  error  of  100%  associated  with  that  of  a  sample 
number  of  unity.  Carrying  this  argument  to  a  200  x  200  grid  of  dimension  4Q0a^  ,  there  will 
be  100  interconnected  pores  corresponding  to  a  statistical  error  of  10%.  Similarly,  a  grid  of 
scale  4000a,nax  will  reduce  the  error  to  1  % . 


Figure  5  illustrates  the  predicted  permeability  measurement  error  associated  with  the 
soil  sample  characterized  above  (a„u„  =0.3  cm).  Note  that  the  errors  associated  with  the 
measurement  of  permeability  become  negligible  as  the  measurement  scale  size  approaches 
several  meters.  This  has  been  confirmed  by  the  infiltration  data  of  Shouse  et.  al.  (1994).  The 
measured  value  of  hydraulic  conductivity  asymptotes  to  2  cm/hr  at  measurement  scales 
greater  than  4  meters.  At  smaller  measurement  scales,  the  inferred  measurement  error  is 
calculated  under  the  assumption  that  the  asymptote  is  precisely  2  cm/hr.  The  excellent 
agreement  between  the  predicted  and  inferred  error  supports  the  extension  of  the  pore  tree 
model  to  describe  a  porous  permeable  media. 


SCALE  OF  MEASUREMENT  (m) 


Figure  5.  Error  in  Permeability  Associated  with  the  Measurement  Scale 


V.  BULK  GASEOUS  DIFFUSION  IN  A  PARTIALLY  SATURATED  MEDIA 


The  extended  pore  tree  model  is  readily  adapted  to  a  partially  saturated  media  through 
the  assumption  that  all  of  the  water  is  contained  in  pore  sizes  between  r^iio  and  r^at  while  only 
gas  is  contained  between  r^^t  and  r„,ax-  Since,  by  Eq.  (6),  porosity  is  distributed  as  In(rp)  in 
pore  space,  the  porosity  associated  with  the  air  filled  pores  (6J  is  approximated  by 


(33) 


where  r^at  is  treated  as  an  independent  variable  of  the  saturated  pore  structure.  No  gaseous 
diffusion  is  allowed  within  rp  <  rgat-  With  this  restriction,  the  extended  pore  tree  model  is  used 
to  develop  an  explicit  relationship  between  bulk  gaseous  diffusivity  and  the  permeability 
which  is  validated  through  the  diffusivity  data  of  Washington  et.  al.,  (1994).  It  is 
demonstrated  that  the  gas  diffusivity  scales  as  l/r^at  and  it  is  the  sensitivity  of  to  the 
saturated  volume  that  controls  the  saturated  diffusivity. 

The  diffusive  mass  flux  in  a  single  pore  is  given  by 


M(r)  =Jt  D„r- 
P'  8  P 


'-ap  ^ 


ctr 


(34) 


where  Dg  is  the  continuum  gas  diffusion  coefficient  (Dg=0.2  cmVs).  Continuum  gas  phase 
diffusion  is  valid  only  for  pore  radii  sufficiently  large  that  Dg  is  greater  than  the  Knudsen 
diffusion  coefficient,  DK„=2Vrp/3,  where  V  is  the  mean  thermal  speed  of  a  molecule. 
Knudsen  diffusion  is  characterized  by  gas  collisions  with  the  pore  walls  and  is  valid  only  for 

Tp  <  3Dg/2V  =  0.1pm.  The  limit  of  validity  of  Eq.  (34)  is  denoted  by  ,  the  larger  of 
either  r^^,  or  3Dg/2V.  This  limit  restricts  continuum  gaseous  diffusion  from  both  the  saturated 
pores  and  the  unsaturated  pores  controlled  by  free  molecule  flow. 


Just  as  in  the  case  of  convection,  it  must  be  determined  whether  the  mass  flux  across 
plane  AA  in  Fig.  4  is  determined  by  the  interconnectivity  of  the  smaller  pores  out  of  the 
plane  (Fig.  4b)  or  by  only  those  pores  that  are  interconnected  in  the  plane  (Fig.  4a).  Consider 
any  pore  of  radius  r,  in  plane  AA  of  Fig.  4b  to  be  the  trunk  of  a  tree.  Each  pore  of  size  rp 
within  the  tree  possesses  the  probability  of  being  interconnected  and  each 

interconnected  pore  in  the  tree  will  carry  the  mass  flow  rate  Afp(rp)  .  Since  there  are  N,drp 

(Eq.  10:  pores  in  size  range  drp  within  the  tree,  the  total  mass  flow  rate 

through  each  and  every  trunk  of  radius  r^  in  plane  AA  becomes 


dr^ 


(35) 
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Integration  of  Eq.  (35)  yields  the  mass  flux  (case  4b)  for  each  trunk  of  radius  rj 


7cZ)  6  r: 


4r 


sat 


(36) 


where  a  ln(r„ax/rsat)  form  was  eliminated  via  Eq.  (33). 

If  the  mass  flux  through  plane  AA  is  limited  by  the  pores  that  are  interconnected  in 

that  plane  (case  4a),  the  mass  flux  is  expressed  as  PjirJ,  and  it  is  immediately  seen 
that 


(37) 


i.e. ,  case  4b  dominates  case  4a.  Since  the  mass  diffusion  in  plane  AA  is  determined  by  the 
pore  interconnectivity  of  the  smaller  pores  outside  of  plane  AA,  the  saturation  of  those 
smaller  pores  becomes  an  important  element  in  the  bulk  gaseous  diffusion. 

Since  each  pore  of  radius  r,  in  plane  AA  carries  mass  flux  ,  the  bulk  diffusion 

coefficient  is  obtained  by  integrating  Eq.  (36)  over  all  pores  in  that  plane 


or,  upon  integration 


D 


bulk 


I 


4r* 


D  00  r 

n  -  ^  a'max 

^bulk 


sat 


(38) 


(39) 


and  upon  eliminating  r,„^  via  Eq.  (31),  the  bulk  diffusivity  Dtuik  is  expressed  in  terms  of  the 
permeability  k. 


D 


2D 


bulk 


'  sat 


(40) 


The  bulk  diffusivity  cannot  increase  indefinitely  with  increasing  permeability  as 
inferred  by  Eq.  (40).  In  deriving  this  expression,  the  mass  flux  through  the  interconnected 
branches  of  the  tree  was  not  constrained  from  exceeding  the  diffusive  capabilities  of  the  trunk 
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itself.  To  correct  this  potential  problem,  the  limit  of  Dbuik  is  determined  as  the  maximum 
diffusive  flux  in  pore  r,  integrated  over  all  pores  in  plane  AA. 


(41) 

This  limit  is  illustrated  in  Figure  6  together  with  the  predicted  values  of  Dbun-  for  an 
extended  range  of  values  of  permeability  and  the  saturation  radius,  r,at.  Model  predictions 
correspond  to  the  measured  values  of  Q^=0.2  and  6=0.5  from  Washington  et.  al.,  (1994), 
and  the  diffusivity  data  suggest  a  value  of  in  the  range  of  10  pm  to  100  pm.  An  exact 
comparison  of  the  present  theory  to  the  least  squared  fit  of  the  data  suggests  a  value  of 
30  pm.  While  the  excellent  agreement  with  the  data  of  Washington  et.  al.,  (1994)  does 
substantiate  the  present  theory,  there  is  clearly  a  very  wide  range  of  possible  values  for  Dbuik 
in  partially  saturated  soil  which  will  depend  upon  an  unknown  saturation  radius.  A  two  order 
of  magnitude  decrease  in  the  saturation  radius  will  increase  the  bulk  diffusivity  by  two  orders 
of  magnitude  and  yet  the  corresponding  increase  in  the  air  filled  porosity  is,  by  (33),  only 
33%.  Hence,  field  measurements  of  the  unsaturated  volume  are  not  sufficiently  accurate  to 
correlate  bulk  diffusivities.  If  bulk  diffiisivities  are  to  be  correlated  with  field  data,  such 
measurements  should  attempt  to  measure  the  saturation  radius. 


10”  10"  lO"  10^  10^  10'  10"*  10^  10"* 


PERMEABILITY  (cm') 


Figure  6.  Bulk  Diffusivity  of  a  Partially  Saturated  Soil 


VI.  PERMEABLE  SUB-RANGE 


The  present  theory  offers  the  opportunity  to  investigate  some  of  the  features  of  the 
pore  structure  that  are  relevant  to  the  fine  scale  contaminant  transport  and  remediation 
reactions.  One  such  feature  is  the  distribution  of  permeability  in  pore  space.  While  the  bulk 
permeability  is  dominated  by  the  largest  pores,  some  permeability  occurs  in  smaller  pores.  It 
is  the  balance  of  the  sub-scale  convection  with  the  small  scale  diffusion  that  will  control 
contaminant  transport  and  in-situ  remediation. 

The  analysis  and  data  comparison  on  the  errors  introduced  by  the  measurement  scale 
size,  section  IV,  has  shown  that  soil  with  a  bulk  permeability  of  0.6  Darcy  possesses,  on 
average,  one  large  interconnected  pore  in  any  12  cm  by  12  cm  cross  section.  To  allow  the 
penetration  of  the  permeate  throughout  such  a  coarse  grid,  it  will  be  shown  that  there  is  a 
very  extensive  permeable  sub-range  in  the  pore  structure  which  does  not  contribute 
significantly  to  the  bulk  permeability  but  in  which  convection  dominates  diffusion. 

To  evaluate  the  permeable  sub-range,  Eq.  (31)  is  rewritten  to  express  the  permeability 
contained  in  pores  of  radius  rp  (rp  <  r^^^x)  or  smaller 


K  = 


16  P 


^2 


(42) 


where  Iq,  represents  the  permeability  in  the  sub-range  dimension  to  be  defined  by  L„(rp).  It 
was  deduced  above  that  L^/fp)  must  be  selected  such  that  there  is  at  least  one  interconnected 
pore  of  radius  Tp  in  cross-sectional  area  'Lj  in  order  for  k^^  to  be  accurate  to  order  unity. 
Subsequently,  a  20  x  20  grid  of  spheres  of  radius  a„  will  contain  one  interconnected  pore  of 
radius  rp 


L  (r  )  =40a 


where  a^^  is  given  by  Eq.  (32)  with  k  replaced  with  k,,. 


(43) 


(44) 


Eqs.  (42)  to  (44)  represent  the  scaling  of  the  pore  structure  with  scale  size  for  all  rp  <  r„,gx  or 
for  all  a  <  a^^,;  and  all  L^/fp)  ^  L„(rmax)  •  The  balance  of  convection  and  diffusion  within 
these  scale  sizes  may  be  assessed. 

The  convective  flux  of  species  "c"  with  mass  density  p,,  across  area  Lj  is  denoted  by 


(45) 
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and  the  diffusive  mass  flux  across  the  same  cross  section  is  expressed  as 


H6) 

where  the  characteristic  gradient  of  p^.  is  pj  L„  and  the  effective  diffusion  coefficient  is 

The  balance  between  convection  and  diffusion  occurs  on  the  length  scale  L  where  the 
convective  and  diffusive  mass  fluxes  are  equal.  From  Eqs.  (45)  and  (46) 


L  =  100 


\l/3 


[e^(-dpjdx) 


(47) 


The  boundary  between  convective  and  diffusive  transport  is  illustrated  in  Figure  7  for  a  solute 
in  a  saturated  porous  media.  The  effective  diffusion  coefficient  is  assumed  to  be  given  by  the 
porosity  times  the  solute  self  diffusion  coefficient  D,  (D,=  10'^  cm^/s).  The  boundary  between 


Figure  7.  Permeable  Sub-Range  for  Convection  in  a  Saturated  Porous  Media 
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sub-scale  convection  and  diffusion  is  illustrated  for  the  solute  in  water  (|i  =  0.01  poise)  at 
pressure  gradients  of  0.01  atm/m,  0.1  atm/m  (hydraulic)  and  10  atm/m.  Also  illustrated  is  the 
length  scale  associated  with  the  spacing  of  the  large  pores  controlling  the  bulk  permeability 
(Eqs.  (43)  and  (44)  in  the  limit  of  ->  ).  The  results  depict  a  very  extensive  permeable 

sub-range  at  bulk  permeabilities  characteristic  of  soil,  sand  and  gravel.  For  soil  with  a  bulk 
permeability  of  1  Darcy,  the  large  interconnected  pores  are  nominally  15  cm  apart.  The  soil 
matrix  between  the  large  pores  experiences  a  convective  mass  flux  through  its  smaller  pores 
whose  mass  flow  rates  are  much  smaller  than  that  of  the  bulk  permeability  but  greater  than 
that  possible  by  diffusion.  At  length  scales  of  order  one  centimeter,  diffusion  becomes  rate 
limiting.  Hence,  this  pore  structure  model  provides  a  methodology  for  determining  the  length 
scale  separating  the  mobile  and  immobile  regions  of  the  soil.  The  theory  depicts  that  the  size 
of  the  immobile  region  is  primarily  a  function  of  the  pressure  gradient. 

When  the  above  analysis  is  applied  to  gas  phase  diffusion  in  an  unsaturated  media, 
p.  is  of  the  order  of  10  "^  p>oise  and  is  0.05  cm^/s  (Figure  6).  The  boundary  between 
convection  and  diffusion  occurs  on  length  scales  five  times  as  large  as  those  depicted  in 
Figure  7  for  the  same  pressure  gradients.  While  the  size  of  the  permeable  sub-range  is 
severely  reduced,  the  entire  distribution  of  species  "c"  is  still  limited  by  a  diffusive  length 
scale  that  is  primarily  a  function  of  the  gas  pressure  gradient.  Hence,  the  length  scale 
separating  the  mobile  and  immobile  regions  of  the  soil  is  of  the  order  of  five  centimeters  at 
gas  pressure  gradients  equal  to  the  hydraulic  pressure  gradient  (0.1  atm/m)  and  of  the  order 
of  10  cm  at  0.01  atm/m. 

Similarly,  when  the  above  analysis  is  applied  to  gas  phase  diffusion  in  a  partially 
saturated  media,  p  is  of  order  lO  '^  poise  and  may  be  as  low  as  10  *  cmVs  (Figure  6). 
Under  these  conditions,  the  boundary  between  convection  and  diffusion  occurs  at  length  scales 
identical  to  those  depicted  in  Figure  7  and  the  existence  of  a  permeable  sub-range  is 
preserved.  However,  for  a  partially  saturated  media,  Db„|,,  should  be  used  for  .  From  Eq. 
(40)  we  write 


rsat 


(48) 


from  which  it  follows  that  the  length  scale  separating  the  mobile  and  immobile  regions  of  a 
partially  saturated  soil  has  a  slightly  stronger  pressure  gradient  dependence 


(  p  6  ^  D  0 

L=44\  o  s  ^  (49) 

r^(-dpldx)^ 

It  is  anticipated  that  the  concept  of  the  permeable  sub-range  will  help  researchers 
develop  relatively  simple,  physics  based  submodels  for  the  ground  water/remediation  codes.  If 
these  submodels  were  tested  independently  from  the  codes,  the  codes  themselves  would  not 
require  parameter  "fitting"  and  could  become  more  directive  than  interactive. 
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VII.  SUMMARY 


The  pore  tree  model  has  been  extended  to  describe  the  permeable  pore  structure  which 
characterizes  the  subsurface  transport  of  gas  and  water  in  soil,  the  dispersion  of  contaminants, 
and  the  in-situ  remediation  of  contaminated  sites.  The  random  nature  of  the  pore  structure, 
which  formed  the  basis  of  the  statistical  derivation  of  the  pore  tree,  is  applied  to  porous  soil 
and  sand.  The  interconnectivity  of  the  pore  structure  is  obtained  via  a  statistical  determination 
of  the  "branches"  that  are  common  to  several  trees  to  allow  convection  and  bulk  diffusion 
through  the  large  scale  (mobile)  structure  in  addition  to  diffusion  and  coupled  chemical 
reactions  within  the  smaller  scale  (immobile)  structure.  The  statistical  analysis  reported  above 
has  determined  that  the  probability  of  pore  interconnectivity  extends  across  the  entire  pore 
size  range,  with  a  slight  increase  in  the  probability  accompanying  a  decreasing  pore  size. 
While  permeability  is  dominated  by  the  largest  pores,  it  is  also  important  to  establish  the  level 
of  convection  and  diffusion  that  is  occurring  at  the  intermediate  scales  in  order  to  accurately 
relate  large  scale  bulk  transport,  small  scale  diffusion  and  coupled  chemical  reactions. 

The  permeability  across  a  given  plane  is  limited  by  the  largest  pores  that  are 
interconnected  in  that  plane.  The  statistical  analysis  has  determined  that  approximately  one 
quarter  of  one  percent  of  all  large  pores  are  interconnected.  This  establishes  a  very  coarse 
grid  for  the  permeability  which  leads  to  measurement  scale  size  errors.  The  extended  pore 
tree  model  has  successfully  explained  the  measurement  errors  in  the  permeability  of  soil  due 
to  the  measurement  scale  size  (Shouse,  et.al.,  1994)  which  has  indirectly  confirmed  the  low 
probability  of  the  interconnectivity. 

The  bulk  gaseous  diffusivity  across  a  given  plane  is  shown  to  be  limited  by  the 
interconnectivity  of  the  smaller  branches  outside  of  that  plane.  These  small  pores  may  be 
saturated,  resulting  in  a  strong  dependence  of  the  diffusivity  on  the  radius  of  the  saturated 
pore.  A  comparison  of  the  present  theory  to  the  diffusivity  data  of  Washington  et.  al.,  (1994) 
suggests  a  saturation  radius  of  30  pm.  \^ile  the  excellent  agreement  with  the  data  does 
substantiate  the  present  theory,  the  diffusivity  in  partially  saturated  soil  is  very  sensitive  to  an 
unknown  saturation  radius.  If  bulk  diffusivities  are  to  be  correlated  with  field  data,  such 
measurements  should  attempt  to  measure  the  saturation  radius. 

The  permeability  and  the  bulk  diffusivity  have  tested  two  extreme  limits  of  the  pore 
structure  and  pore  interconnectivity  concepts.  Permeability  is  limited  by  the  in  plane 
interconnectivity  (Fig.4a)  and  bulk  gaseous  diffusion  is  limited  by  the  out  of  plane  (Fig.  4b) 
interconnectivity.  Permeability  is  limited  by  the  large  pore  interconnectivity  and  bulk  diffusion 
is  limited  by  the  interconnectivity  of  the  smaller  pores.  The  apparent  success  of  these 
concepts  over  a  very  broad  pore  size  range  suggests  that  the  extended  pore  tree  model  will 
accurately  describe  the  relationship  between  large  scale  convection  and  small  scale  diffusive 
transport. 

Analysis  of  the  permeability  has  utilized  the  interconnectivity  of  the  pores  to  determine 
the  distribution  of  the  permeability  with  f>ore  size.  This  analysis  suggests  the  existence  of  a 
permeable  sub-range  in  the  pore  structure  which  does  not  contribute  significantly  to  the  bulk 
permeability  but  in  which  convection  dominates  diffusion.  It  is  the  balance  of  the  sub-scale 
convection  with  the  small  scale  diffusion  that  will  control  contaminant  transport  and 
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in-situ  remediation.  Preliminary  estimates  suggest  that  the  length  scale  separating  the  mobile 
and  immobile  regions  of  the  soil  is  of  the  order  of  one  cm.  Smaller  grains  are,  on  a  unit  mass 
or  size  basis,  more  reactive  than  the  larger  grains.  Hence  smaller  grains  will  locally  deplete 
more  nutrients  or  remediation  chemicals  from  the  convective  flow  than  their  larger 
counterparts.  Simultaneously,  the  smaller  pores  surrounding  the  smaller  grains  offer  more 
resistance  to  convection  and  a  reduced  amount  of  nutrients  or  remediation  chemicals  will  be 
available  in  the  convective  flow.  The  size  distribution  of  the  pores  and  grains,  and  the 
variations  in  fluid  velocity  within  and  between  pores  of  different  sizes  is  critical  to  interfacing 
these  transport  processes.  Future  plans  include  the  development  of  a  methodology  to  couple 
convection  and  small  scale  diffusion  upon  which  coupled  chemical  reactions  may  be  added  to 
accurately  describe  contaminant  transport  and  in-situ  remediation. 
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